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Top-quark pair + 1-jet production at next-to-leading order QCD 
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Top-quark pair production with an additional jet is an important signal and background process at the LHC. We present the 
next-to-leading order QCD calculation for this process and show results for integrated as well as differential cross sections. 



1. Introduction 

With a mass of 172.6 GeV the top-quark is by 
far the heaviest elementary fermion in the Stan- 
dard Model (SM). Its mass is more than 30 times 
larger than the mass of the next heaviest fermion, 
the bottom-quark. The large mass has lead to vari- 
ous speculations whether the top-quark behaves as a 
normal quark or whether it plays a special role. In 
particular, the fact that the top-quark mass is close 
to the scale of electroweak symmety breaking — or 
equivalently, that the Yukawa coupling to the Higgs is 
very close to one — has motivated different scenarios 
in which the top-quark drives the electroweak symme- 
try breaking. For a recent overview we refer the inter- 
ested reader to Refs. ||1|2| . In the context of the SM 
the top-quark interactions are completely determined 
through the gauge structure. The only free param- 
eter appearing in top-quark physics is the top-quark 
mass. Once this parameter is measured all remain- 
ing properties are predicted. An important task for 
the ongoing Tevatron collider and the upcoming LHC 
is the precise measurement of the top-quark proper- 
ties. In this context the production of a top-quark pair 
togehter with an additional jet is an important reac- 
tion. This becomes already clear from the simple ob- 
servation that a substantial number of events in the 
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inclusive top-quark sample is accompanied by an ad- 
ditional jet. Depending on the energy of the additional 
jet the fraction of events with an additional jet can eas- 
ily be of the order of 10-30%. For a more precise un- 
derstanding of the topology of top-quark events it is 
thus important to have also an improved understand- 
ing of top-quark pair production together with a jet. 
Recently it has been argued that the inclusive top- 
quark pair cross section can be measured at the LHC 
with an experimental uncertainty of about 5%. The 
available theoretical predictions at next-to-leading or- 
der (NLO) togehter with resummation lead to a the- 
oretical uncertainty of more than 10% II3I4I5I . The 
dominant uncertainy comes from the residual scale 
dependence, while the uncertainties due to the par- 
ton distribution functions are comparably small. In 
Ref. |3 1 parts of the full next-to-next-to-leading order 
(NNLO) were derived from general arguments. It has 
been shown that an NNLO calculation will reduce the 
theoretical uncertainty down to 5%. Apart from the 
two-loop corrections (where progress has been made 
recently B6I7II ) one important ingredient in going to 
NNLO are the one-loop coiTections to tt+ 1-jet pro- 
duction. Apart from its significance as signal process 
it turns out that tt+ 1-jet production is also an im- 
portant background to various new physics searches. 
A prominent example is Higgs production via vector- 
boson fusion. This reaction represents an important 
discovery channel for a SM Higgs boson in the mass 
range of 120-180 GeV. The major background to this 
reaction is due to tt+ 1-jet H), again underlining the 
need for precise theoretical predictions for this pro- 
cess. It is well known that predictions at leading-order 
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(LO) in the coupling constant of QCD are plagued by 
large uncertainties. In many cases the LO predictions 
in QCD give only a rough estimate. Only by includ- 
ing NLO corrections a quantitative reliable prediction 
can be obtained. Given that the conceptual problems 
of doing such calculations are solved since quite some 
time, one might think that doing the required calcula- 
tions should be a straightforward task. Unfortunately 
it turns out that this is not the case. The calcula- 
tion of QCD corrections for 2^3 and 2-^4 reac- 
tions is still a highly non-trivial task — not speaking 
about reactions with an even higher multiplicity. In 
general the problem can be attributed to the fact that 
the corrsponding matrix elements are complex func- 
tions of many variables so that an analytic treatment 
is no longer feasible due to the large size of the ex- 
pressions. A solution to this problem is to resort to 
numerical methods. While in principle fine one is in 
many cases plagued by numerical instabilities and the 
long runtime. In particular, the reduction of one-loop 
tensor integrals to scalar one-loop integrals is in gen- 
eral difficult to do in a numerically stable way. In that 
context the calculation of the one-loop corrections to 
top-quark pair production with an additional jet is also 
interesting as a benchmark process for the develope- 
ment of new methods. In the present article we will 
briefly comment on the calculation of the NLO cor- 
rections. In addition we will show results for inte- 
grated quantities as well as for differential distribu- 
tions. 

2. Calculation 

2.1. Born approximation 

In Born approximation the partonic reactions are 
88 -^ %. qq -^ %. 18 -^ "?. and gq -^ ttq. The 
last three reactions are related by crossing. In Born 
approximation various well-tested approaches to cal- 
culate the required matrix elements exist. We used 
recurrence relation a la Berends and Giele f^l and 
a Feynman-diagram-based approach. In both cases 
four-dimensional helicity methods were employed. In 
addition we also used Madgraph |10| for checking. 
We found complete agreement of the different meth- 
ods. The explicit LO calculation shows indeed the 
large scale dependence as expected. Without going 
into details we just mention that in LO the importance 
of the individual partonic channels does not follow 



the pattern known from inclusive top-quark pair pro- 
duction. While at the Tevatron the situation tt+ 1-jet 
production is similar to the inclusive reaction, that is 
the total cross section is dominated by the quark-anti- 
quark channel followed by the gluon-fusion process, 
the situation at the LHC is different from inclusive 
production. The most important channel is gluon fu- 
sion, but in contrast to inclusive production the sec- 
ond important channel is the q'g-channel. This is due 
to the large parton luminosity for this channel and due 
to a sizeable partonic cross section. 

2.2. Virtual corrections 

The virtual corrections consist of the one-loop cor- 
rections to the LO reactions. One can classify the cor- 
rections into self-energy corrections, vertex correc- 
tions, box-type corrections, and pentagon-type cor- 
rections where all the external legs are connected to 
one loop thus forming a pentagon. The latter are the 
most complicated ones owing to their complexity and 
the involved tensor integrals. The challenging step 
in this context is the numerically fast and stable re- 
duction of the tensor integrals to scalar one-loop in- 
tegrals. To ensure the correctness of our results we 
did two independent calculations of the virtual cor- 
rections using as far as possible different methods and 
also different tools. In one calculation the virtual cor- 
rections are essentially obtained following the proce- 
dure described in Ref. |11|, where tt + H production 
at hadron colliders was considered. Feynman dia- 
grams and amplitudes have been generated with the 
FeynArts package 1 12 13 1 and further processed with 
in-house Mathematica routines, which automatically 
create an output in Fortran. The IR (soft and collinear) 
singularities are analytically separated from the finite 
remainder as described in Refs. II11II4II . The tensor 
integrals appearing in the pentagon diagrams are di- 
rectly reduced to box integrals following Ref. ifTSJI . 
This method does not introduce inverse Gram deter- 
minants in this step, thereby avoiding notorious nu- 
merical instabilities in regions where these determi- 
nants become small. Box and lower-point integrals 
are reduced a la Passarino-Veltman fW\ to scalar inte- 
grals, which are either calculated analytically or using 
the results of Refs. M17II8I19I . Sufficient numerical 
stability is already achieved in this way. Nevertheless 
the integral evaluation is currently further refined by 
employing the more sophisticated methods described 
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in Ref. ||201 in order to numerically stabilize the tensor 
integrals in exceptional phase-space regions. 

In the second calculation the evaluation of loop 
diagrams starts with the generation of diagrams and 
amplitudes via QGRAF 121], which are then further 
manipulated with Form I22| and automatically trans- 
lated into C++ code. The reduction of the 5-point 
tensor integrals to scalar integrals is performed with 
an extension of the method described in Ref. ||23I . 
In this procedure also inverse Gram determinants of 
four four-momenta are avoided. The lower-point ten- 
sor integrals are reduced using an independent imple- 
mentation of the Passarino-Veltman procedure. The 
IR-finite scalar integrals are evaluated using the FF 
package II24I25I . 



2.3. Real corrections 

The matrix elements for the real corrections are 
given by — > ^^gggg, — > ttqqgg, -^ Uqqq' q' and 
-^ ttqqqq. The various partonic processes are ob- 
tained from these matrix elements by all possible 
crossings of light particles into the initial state. 

The evaluation of the real-emission amplitudes is 
again performed in two independent ways. Both 
evaluations employ the dipole subtraction formalism 
Ii26 27 28 1 for the extraction of IR singularities and 
for their combination with the virtual corrections. 

One calculation of the real corrections results from 
a fully automated calculation based on helicity ampli- 
tudes, as described in Ref. 1291 . Individual helicity 
amplitudes are computed with the help of Berends- 
Giele recurrence relations 1 9 1 . The evaluation of color 
factors and the generation of subtraction terms is au- 
tomated. For the channel gg -^ iigg a dedicated soft- 
insertion routine [30] is used for the generation of 
the phase space. The second calculation uses for 
the LO 2^3 processes and the gg -^ iigg process 
optimized code obtained from a Feynman diagram- 
matic approach. As in the calculation described be- 
fore standard techniques like color decomposition and 
the use of helicity amplitudes are employed. For the 
2^4 processes including light quarks, Madgraph 
ifTOl has been used. The subtraction terms accord- 
ing to Ref. 1 28 1 are obtained in a semi-automatized 
manner based on an in-house library written in C++. 



3. Results 

In the following we consistently use the CTEQ6 
113 11321 set of parton distribution functions (PDFs). In 
detail, we take CTEQ6L1 PDFs with a 1-loop run- 
ning tts in LO and CTEQ6M PDFs with a 2-loop 
running tts in NLO. The number of active flavours 
is Np = 5, and the respective QCD parameters are 
A^° = 165 Me V and Ap == 226 MeV. Note that the 
top-quark loop in the gluon self-energy is subtracted 
at zero momentum. In this scheme the running of tts 
is generated solely by the contributions of the light 
quark and gluon loops. The top-quark mass is renor- 
malized in the on-shell scheme, as numerical value 
we take Wt = 174GeV. 

We apply the jet algorithm of Ref. 1331 with R = 
1 for the definition of the tagged hard jet. Unless 
stated otherwise we require a transverse momentum 
of PTjet > P™' = 20GeV for the hardest jet. The out- 
going (anti-)top-quarks are neither affected by the jet 
algorithm nor by the phase-space cut. Note that the 
LO prediction and the virtual corrections are not influ- 
enced by the jet algorithm, but the real corrections are. 

In Fig. [T]the scale dependence of the NLO cross 
sections is shown. For comparison, the LO results 
are included as well. Varying the scale in the usual 
range that is a factor 2 up and down around a cen- 
tral scale we observe a reduction of the scale depen- 
dence of a factor 3 at Tevatron and about a factor 6 at 
the LHC. From the absolute values of the cross sec- 
tions one concludes that indeed a significant contri- 
bution to the inclusive top-quark pair cross section 
comes from tt+ 1-jet. For reasonable scale choices 
one can also observe that the NLO corrections are of 
moderate size. In particular around /j — nit the cor- 
rections are small. The interference of C-odd contri- 
butions of the amplitude with C-even parts, where C 
denotes the charge conjugation, produces a forward- 
backward charge asymmetry at the Tevatron. For the 
inclusive top-quark pair sample, where this asym- 
metry is an NLO effect, this has been studied in 
Refs. I135I36I371 . In the case of tt + 1-jet this asym- 
metry appears already in LO. In LO the asymmetry is 
defined by 
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where y^ denotes the rapidity of the top-quark. Cross- 
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Figure 1. Scale dependence of the LO and NLO cross sections for tt + 1-jet production at the Tevatron (left) and 
the LHC (right) as taken from Ref. [i34J . with the renormalization scale (pr) and the factorization scale (/j/) set to /j. 



section contributions (3{yi < 0) correspond to top- 
quarks in the forward or backward hemispheres, re- 
spectively, where incoming protons fly into the for- 
ward direction by definition. Denoting the corre- 
sponding NLO contributions to the cross sections by 
SOj^LO' ^^ define the asymmetry at NLO by 
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i.e. via a consistent expansion in ttj. Note, however, 
that the LO cross sections in Eq. Q are evaluated in 
the NLO setup (PDFs, tts). The results for the asym- 
metry for different scale choices are shown in Fig.|2] 
At LO we find an asymmetry of about —8%. The 
scale dependence is rather small. This is a conse- 
quence of the fact that ttj cancels exactly between the 
numerator and the denominator In addition the resid- 
ual factorization scale dependence also cancels to a 
large extent in the ratio. At NLO we find a large cor- 
rection compared to the LO result. The asymmetry 
is almost washed out at NLO. The scale dependence 
is increased in NLO which seems natural given the 
small dependence in LO. To investigate the origin of 
the large NLO corrections to the asymmetry we stud- 
ied the dependence on p™*, the minimal p-Y used to 
resolve the additional jet. The results are shown in 
Tab. [U A strong dependence of the cross section on 
p™' is observed. For all p™' values we find that the 
NLO corrections to the cross section are of moderate 
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Figure 2. Scale dependence of the LO and NLO 
forward-backward charge asymmetry of the top- 
quark in pp -^ tt+jet+X at the Tevatron as taken from 
Ref. 1 34 1 with /J = /J/ =/Jr. 
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Table 1 

Cross section and forward-backward charge asymmetry at the Tevatron for different values of p™' used to define 
the minimal transverse momentum pj of the additional jet (jj = /jf = /Jr = nii). The upper and lower indices are 
the shifts towards /j = m,/2 and /j = 2mi. 





cross section [pb] 


chai-ge asymmetry [%] 


pT [GeV] 


LO 


NLO 


LO 


NLO 


20 


1.583(2)1JJ:^^ 


i.79i(i)i!J:^t 


-7.69(4)l«:i«5 


-1.77(5)10:^3« 


30 


0.984(1)1°^° 


1.1194(8)1!^:^^ 


-8.29(5)l«:i25 


-2.27(4)1°:^! 


40 


0.6632(8)1^:^1 


0.7504(5)+!5?f 


-8.72(5)1«:|^ 


-2.73(4)1!J:^5 


50 


0.4670(6)155:2^ 


0.5244(4)1«:«^^ 


-8.96(5)l«|t 


-3.05(4)1«« 



size. While the cross section evidently has a large sen- 
sitivity on p™', the dependence of the asymmetry on 
/7™' is slightly less pronounced. On the other hand we 
find, independent of the chosen p™', a significant dif- 
ference between LO and NLO for the asymmetry in- 
dicating that the definition of the asymmetry Eq. (|2]i is 
not stable with respect to higher-order corrections in- 
dependent of the p™' value. In Fig.[3]the rapidity dis- 
tribution for the top-quark is shown. One can easily 
observe the asymmetry of the LO distribution leading 
to the aforementioned —8% for the integrated asym- 
metry. The NLO corrections are different in shape 
compared to LO. They are larger in the forward di- 
rection than in the backward direction. The result is 
that the asymmetry is washed out by the NLO correc- 
tion. In the lower plot of Fig. [3] the /T-factor and the 
scale dependence is shown. The corrections are be- 
tween +20 and +40%. The large scale dependence 
of the LO result is again significantly improved when 
the NLO corrections are taken into account. In Fig.|4] 
the px-distribution of the top-quark is shown — now 
for the LHC. For not too large pj the corrections are 
of moderate size. Starting with +30 % for small pj 
and reaching —50 % at pj « 450 GeV, indicating that 
the /T-factor is strongly phase-space dependent. Up to 
a Pj of about 450 GeVwe also find an important im- 
provement of the scale uncertainty. 

Conclusions: Predictions for tt+jet production at 
hadron colliders have been reviewed at NLO QCD. 
For the cross section the NLO corrections drastically 
reduce the scale dependence of the LO predictions, 
which is of the order of 100%. The charge asym- 
metry of the top-quarks, which is measured at the 
Tevatron, is significantly decreased at NLO and is al- 
most washed out by the residual scale dependence. In 
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Figure 3. Rapidity distribution of the top-quark. In 
the upper plot the thick (thin) line shows the NLO 
(LO) result. In the lower plot the thick line shows the 
/T-factor, the two bands show the scale dependence 
in LO (outer band) and NLO (inner band) when the 
scale (fi — pr— Pf) is varied by factor 2 up and down 
around the central value p — mi. 
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Figure 4. Same as Fig.|3]for the px -distribution of the 
top-quark at LHC. 



addition we have studied the /:i™'-dependence of the 
asymmetry in NLO. Further refinements of the pre- 
cise definition of the charge asymmetry are required 
to stabilize the predictions with respect to higher- 
order corrections. First resuhs for differential distri- 
butions have been presented. The corrections are well 
under control over a large phase-space region and the 
scale uncertainty is again improved compared to the 
LO results. 
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